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ABSTRACT 

Accurate determinations of stellar mass functions and ages of stellar populations are crucial to much of 
astrophysics. We analyse the evolution of stellar mass functions of coeval main sequence stars including all 
relevant aspects of single- and binary-star evolution. We show that the slope of the upper part of the mass 
function in a stellar cluster can be quite different to the slope of the initial mass function. Wind mass loss 
from massive stars leads to an accumulation of stars which is visible as a peak at the high mass end of mass 
functions, thereby flattening the mass function slope. Mass accretion and mergers in close binary systems 
create a tail of rejuvenated binary products. These blue straggler stars extend the single star mass function by 
up to a factor of two in mass and can appear up to ten times younger than their parent stellar cluster. Cluster 
ages derived from their most massive stars that are close to the turn-off may thus be signiflcantly biased. To 
overcome such difficulties, we propose the use of the binary tail of stellar mass functions as an unambiguous 
clock to derive the cluster age because the location of the onset of the binary tail identifies the cluster turn-off 
mass. It is indicated by a pronounced jump in the mass function of old stellar populations and by the wind mass 
loss peak in young stellar populations. We further characterise the binary induced blue straggler population in 
star clusters in terms of their frequency, binary fraction and apparent age. 

Subject headings: (stars:) binaries: general — (stars:) blue stragglers — stars: luminosity function, mass 
function — stars: mass-loss 


1. INTRODUCTION 

Stellar mass functions are important for population studies, 
both nearby and at high redshift. Mass functions are not static 
but change their shape when the frequencies of stars or their 
masses are altered over time (Scalo 1986; Kroupa et al. 2013). 
Possible causes include the evaporation of stars and mass seg¬ 
regation in star clusters (de Grijs et al. 2002; McLaughlin & 
Fall 2008; Harfst et al. 2010; Habibi et al. 2013) and mass loss 
in the course of stellar evolution (Danger 2012). These mecha¬ 
nisms leave characteristic fingerprints in mass functions from 
which insights into the evolutionary status of stellar popula¬ 
tions and the mechanisms themselves can be gained. 

Massive stars are subject to strong stellar wind mass loss, 
which decreases their masses already on the main sequence 
(MS), directly affecting mass functions. Furthermore, it 
emerges that binary stars play an important role in stellar pop¬ 
ulations of various ages and even dominate the evolution of 
massive stars (Sana et al. 2012). The multiplicity fraction, i.e. 
the number of multiple stars divided by the total number of 
stellar systems, is larger for higher masses (e.g. Bastian et al. 
2010; Duchene & Kraus 2013): it exceeds 40% for solar like 
stars (F, G and K stars) and 70% for the most massive stars 
(0-stars). In close binaries, mass is exchanged between the 
binary components during Roche lobe overflow (RLOF) or in 
stellar mergers, directly affecting stellar masses and hence the 
mass function. The mass gainers are rejuvenated and can ap¬ 
pear much younger than they really are. Some mass gainers 

fabian. Schneider @ physics . ox. ac .uk 

^ Argelander-Institut fur Astronomic der Universitat Bonn, Auf dem 
Hugel 71, 53I2I Bonn, Germany 

^ Department of Physics, University of Oxford, Denys Wilkinson Build¬ 
ing, Keble Road, Oxford 0X1 3RH, UK 

^ Institute of Astronomy, University of Cambridge, Madingley Road, 
Cambridge CBS OHA, UK 

"^Astronomical Institute ’Anton Pannekoek’, Amsterdam University, 
Science Park 904, 1098 XH, Amsterdam, The Netherlands 


may also be visible as blue straggler stars (Braun & Danger 
1995; Dray & Tout 2007; Geller & Mathieu 2011). 

The determination of stellar ages is a fundamental task in 
stellar astrophysics (Soderblom 2010) that can be biased by 
rejuvenated binary products. Various methods are used to de¬ 
termine stellar ages. The surface properties of individual stars 
can be compared with evolutionary model predictions, e.g., 
luminosity, surface temperature or surface gravity (Holmberg 
et al. 2007; Schneider et al. 2014a), rotation rates in low mass 
stars (Barnes 2007) or surface nitrogen abundances in massive 
stars (Kohler et al. 2012). In star clusters, the most widely 
used age determination method compares the main sequence 
turn-off with theoretical isochrones (e.g. Naylor & Jeffries 
2006; Monteiro et al. 2010). Close binary evolution such as 
mass transfer and stellar mergers leads to spurious or inaccu¬ 
rate results in all of these methods. To derive unambiguous 
age estimates one must be able to distinguish between rejuve¬ 
nated binary products and genuinely single stars. 

In this paper we investigate how the modulation of stellar 
mass functions by single and binary star evolution can be used 
both to identify binary products and to derive unambiguous 
stellar ages. More generally, we explore what can be learned 
about stellar evolution and the evolutionary status of whole 
stellar populations from observed mass functions. We also 
investigate quantitatively how single and binary star evolution 
influence the determination of initial mass functions. To that 
end, we perform detailed population synthesis calculations of 
coeval stellar populations using a rapid binary evolution code. 

We describe our method in Sec. 2 and present mass func¬ 
tions of coeval single and binary star populations of ages rang¬ 
ing from 3Myr to IGyr in Sec. 3. Through binary interac¬ 
tions, blue straggler stars are formed in our models. We char¬ 
acterise their binary fraction and ages, and compare their pre¬ 
dicted frequencies to those found in Galactic open star clus¬ 
ters in Sec. 4. Blue straggler stars predominantly populate the 
high mass end of mass functions and may bias determinations 
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of stellar cluster ages. We show how to use mass functions to 
overcome such biases when determining cluster ages in Sec. 5 
and conclude in Sec. 6. 

2. METHOD 

We compute the evolution of single and close, interacting 
binary stars and construct, at predefined ages, mass functions 
by counting how many stars of certain masses exist. This ap¬ 
proach ensures that we factor in all the relevant single and bi¬ 
nary star physics to investigate their infiuence on the present- 
day mass function (PDMF). 

The initial parameter space of binaries is large compared 
to that of single stars, which essentially consists of the initial 
mass. In binaries, there are two masses, the orbital separa¬ 
tion, the eccentricity of the orbit and the relative orientation 
of the spin axis of both stars. We apply some standard sim¬ 
plifications to reduce this huge parameter space: we impose 
circular orbits and that both stellar spins are aligned with the 
orbital angular momentum. All our models are calculated for 
a metallicity Z = 0.02 (unless stated otherwise). Furthermore 
we focus on the main sequence because stars spend typically 
90% of their life in this evolutionary stage. Still, we need to 
follow the evolution of a large number of stellar systems to 
sample the remaining binary parameter space and to resolve 
effects at the high mass end of the PDMF. Hence, we work 
with a rapid binary evolution code. 

2.1. Rapid binary evolution code 

We use the binary population and nucleosynthesis code of 
Izzard et al. (2004, 2006, 2009) with modifications due to de 
Mink et al. (2013) which is based on the rapid binary evolu¬ 
tion code of Hurley et al. (2002). This code uses analytic for¬ 
mulae (Hurley et al. 2000) fitted to detailed single star evolu¬ 
tionary models (Pols et al. 1998) to approximate the evolution 
of single stars for a wide range of masses and metallicities. 

The fitting formulae of Hurley et al. (2000) are based on de¬ 
tailed stellar model sequences of stars with mass up to 50 Mq 
(P ols et al. 1998). The evolution of stars with mass in excess 
of 50 Mq is thus based on extrapolations of the original fitting 
formulae. The MS lifetime, tms, is particularly inaccurately 
extrapolated by the appropriate fit, so we replace it with a 
logarithmic tabular interpolation of the MS lifetimes taken di¬ 
rectly from the models of Pols et al. (1998) in the mass range 
20 < M < 5OM0. More massive than this we extrapolate the 
final two masses in the grid of detailed models (Pols et al. 
1998). This results in a reduction of the MS lifetime of e.g. 
a 100 Mq star at Z = 0.02 from 3.5 to 2.9Myr, i.e. in a re¬ 
duction of 17%, which is in agreement with state-of-the-art 
non-rotating detailed stellar models of Brott et al. (2011) and 
Ekstrom et al. (2012). 

Stellar wind mass loss for stars with luminosities L > 
4000 Lq is given by Nieuwenhuijzen & de Jager (1990). This 
recipe is modified by a factor Z^-^ according to Kudritzki et al. 
(1989) to mimic the impact of the metallicity Z on wind mass 
loss rates. 

Binary stars can exchange mass either by RLOF, wind mass 
transfer or merging. RLOF occurs when one star (hereafter 
the donor) fills its Roche lobe and transfers mass to its com¬ 
panion (hereafter the accretor) through the inner Lagrangian 
point. Depending on the physical state of the donor star one 
distinguishes between Case A, B and C mass transfer. Fol¬ 
lowing the definition of Kippenhahn & Weigert (1967), Case 
A mass transfer occurs during core hydrogen burning. Case B 


after the end of core hydrogen burning and Case C, defined 
by Lauterbom (1970), after the end of core helium burning. 

Our binary evolution code differs from Hurley et al. (2002) 
in its treatment of RLOF. For stable mass transfer it is ex¬ 
pected that the stellar radius R adjusts itself to the Roche lobe 
radius R\^/\.t. R ^ R\^. Whenever RLOF occurs (R > Rl) we 
remove as much mass as needed to shrink the donor star back 
into its Roche lobe. The resulting mass transfer and mass ac¬ 
cretion rates are capped by the thermal timescales of the donor 
and accretor, respectively. 

We follow Hurley et al. (2002) to determine the occurrence 
of common envelope evolution and contact phases that lead to 
stellar mergers. MS mergers are expected either if the initial 
orbital separation is so small that both stars fill their Roche 
lobes and thus come into physical contact or if the mass ratio 
of the accretor to the donor star falls below a certain limit at 
the onset of mass transfer, q = M 2 /Mi < which drives 
the accretor out of thermal equilibrium and hence results in 
a contact system (e.g. Ulrich & Burger 1976; Kippenhahn 
& Meyer-Hofmeister 1977; Neo et al. 1977; Wellstein et al. 

2001) . The critical mass ratio is approximately ^crit,MS = 0.56 
for MS stars (de Mink et al. 2007), ^crit,HG = 0-25 if the donor 
star is a Hertzsprung gap star and is given by a fitting formula 
if the donor star has a deep convective envelope (Hurley et al. 

2002 ) . 

Mass transfer because of either stable RLOF or during a 
stellar merger makes the mass gainers appear younger than 
they really are (Braun & Langer 1995; van Bever & Vanbev- 
eren 1998; Dray & Tout 2007). Such rejuvenated stars may 
stand out as blue stragglers in Hertzsprung-Russell (HR) di¬ 
agrams. Rejuvenation is handled following Tout et al. (1997) 
and Hurley et al. (2002) but with improvements as described 
in Glebbeek & Pols (2008) and de Mink et al. (2013). The 
apparent age T of a star is given by the amount of burnt 
fuel compared to the total available. For a star with MS 
lifetime tms we therefore have (in a linear approximation) 
T = /bumtTMS where /burnt ~ ^^burnt/^^available is the maSS ratio 
of the burnt to totally available fuel. After mass transfer onto 
a MS star with a convective core, the mass of the already burnt 
fuel is given by the fraction of burnt material, /burnt, times the 
convective core mass before mass transfer, M^ (primes indi¬ 
cate quantities after mass transfer). After mass transfer, the 
convective core and hence the available fuel of the accretor 
grow in mass, i.e. M' > Me, because the total stellar mass in¬ 
creases (and vice versa for the donor star). Thus, the fraction 
of burnt fuel after mass transfer is = /bumtMc/M' and the 
apparent age, T' = is 

'Y tit 

~ M' TMS ■ 

This equation holds for the accretor and also for the donor 
when setting Me = M' (no burnt fuel is mixed out of the core 
upon mass loss from the stellar surface) and shows that the 
accreting stars rejuvenate upon mass transfer (T' < T because 
M' > Me and < tms) while the donor stars age (T' > T 
because > tms)- The accretor and donor will have burnt 
less (more) fuel than a single star of the same mass that did 
not accrete (lose) mass. If stars do not have convective cores, 
e.g., stars with initial masses in the range 0.3 < M/Mq <1.3 
or Hertzsprung gap stars which have radiative cores, we set 
Me = M' because no fresh fuel is expected to be mixed into 
their cores. 

To model the rejuvenation of MS mergers we follow de 
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Mink et al. (2013): first, we assume that a fraction /loss of the 
total mass M 3 = Mi + M 2 is lost during the merger; we adopt 
/loss = 0.1. Second, we approximate the core mass fraction /c 
of MS stars according to fitting functions (Glebbeek & Pols 
2008) and estimate the apparent age of the newly formed 
merged star from. 


T3 = tms,3 • 


/c,l 


^MS,l ’ ^MS,2 

Jc,3 


( 2 ) 


where 7i and T 2 are the effective ages, /c,i and /c ,2 are the core 
mass fractions of the progenitor stars and tms denotes the MS 
lifetime of the corresponding star. The denominator contains 
the effective core mass fraction /ff = /c ,3 +/mix * (1 -/c, 3 ) of the 
merger product which is given by its core mass fraction mod¬ 
ified by an additional mixing of /mix = 10 % of the hydrogen- 
rich envelope. The resulting rejuvenation is less than in the 
original Hurley et al. (2002) code, where it is assumed that 
the whole star is mixed, and closer to that seen in hydrody¬ 
namic and smoothed particle hydrodynamics simulations of 
stellar mergers (e.g. Lombardi et al. 1995; Sills et al. 1997, 
2001; Glebbeek & Pols 2008; Glebbeek et al. 2013). 


2 . 2 . Initial distribution functions 

We set up a grid of stellar systems to cover the parameter 
space of single and binary stars and assign each stellar system 
j a probability of existence 5pj. The probabilities 5pj are cal¬ 
culated from 5pj = ^5\nV where is a distribution function 
of the initial masses and the initial orbital periods and ^InV is 
the volume of the parameter space filled by the stellar system 
/ The initial distribution function, reads, 

^ ^ f V^(lnm) single stars, 

\'/(lnmi)/(lnm 2 )x(lnP) binary stars. 


where mi and m 2 are the initial masses of the primary and 
secondary star in binaries, respectively, and P is the initial or¬ 
bital period. The functions '/(Inmi), /(Inm 2 ) and x(lnP) are 
the IMFs of the primary and secondary star and the distribu¬ 
tion function of the initial orbital period, respectively. Stellar 
masses m are given in solar masses. 

Single stars and also the primaries of binary stars, i.e. the 
initially more massive component, are distributed according 
to a Salpeter IMF (Salpeter 1955). The IMF /(InmO is then 

&p &p Y 

'/(Inmi) = —- =mi -—=Ami (4) 

din mi dmi 


with r = -1.35 the slope of the mass function and A the nor¬ 
malization constant such that 



dp 

dmi 


dmi = 1. 


( 5 ) 


The lower and upper mass limits are chosen such that we do 
not exceed the validity of the fitting functions used in the code, 
hence m/ = 0.1 and mu = 100 . 

Duchene & Kraus (2013) review stellar multiplicity (mul¬ 
tiplicity fractions, mass ratio distributions and orbital separa¬ 
tions distributions) and its dependence on primary mass and 
environment. A complete picture is yet lacking (e.g., for MS 
primary stars in the mass range 8-16 Mq) but it seems that a 
fiat mass ratio distribution. 


/(In m 2 ) = 


dp 

din m 2 



oc^, 


( 6 ) 


is reasonable maybe except at the lowest primary masses 
(< IMq). We adopt this as the distribution function of the 
initial mass ratios, meaning that all mass ratios q are equally 
probable (d;?/d^ = const.). We adopt a minimum mass ratio, 

^min — 0.1 /mi. 

In terms of orbital separations, we are only interested in 
close binaries, i.e. binaries that can interact by mass exchange 
at some point during their life. Our binary systems therefore 
have initial orbital separations a between 3 Rq and 10 "^ Rq 
(^ 46 AU). In practice it turns out that only binaries with ini¬ 
tial orbital separations less than about 3 x 10^ Rq ^ 15AU 
interact. Initially wider binaries are effectively single stars. 
The distribution function x(lii^) of the initial orbital periods 
P is given by Sana et al. (2012) for binaries with 0-type com¬ 
panions (i.e. mi > 15) and mass ratios ^ > 0.1 and by a fiat 
distribution in InP for all other binaries, i.e. /(P)dP oc dP/P 
(Opik 1924), 


X(lnP) oc 


const. 


0.15 <logP/d< 3.5 
otherwise. 


( 7 ) 


The lower boundary of the initial orbital separations, ai = 
3Rq, is increased if a star fills its Roche lobe on the zero- 
age main sequence (ZAMS) such that stars cannot interact 
immediately by RLOF. 

For our mass function calculations, we evolve a total of 
2,500,000 binary and 250,000 single stars. The overall 
mass range of 1-100 Mq is subdivided in 10 equally spaced 
mass intervals. Each of the 10 mass intervals covers 2500 
single stars of different initial masses and 100 x 50 x 50 
(mi X qx a) binaries. Masses and orbital separations are dis¬ 
tributed equidistantly on a logarithmic grid while the mass 
ratios are distributed equidistantly on a linear grid. 


2.3. Binary parameter space 

To understand the infiuence of binary interactions on the 
PDMF, we need to know quantitatively how much mass is 
transferred and accreted during binary mass exchange. Given 
an initial primary mass mi, we investigate the binary param¬ 
eter space spanned by the initial secondary mass m 2 {< mi) 
and initial orbital separation a to calculate the amount of 
transferred and accreted mass because of RLOF (AMtrans and 
AMacc respectively) during the MS evolution of the secondary 
star. We further compute the mass transfer efficiency. 


^^acc 


( 8 ) 


during the different mass transfer Cases A, B and C. AMacc 
and AMtrans are then the accreted and the transferred mass, 
respectively, during the different mass transfer cases. 

In Fig. 1 we show the mass transfer efficiency / (Eq. 8 ) in 
binaries with IOMq primary stars as a function of initial mass 
ratio, ^ini, and initial orbital separation, aini. Binaries interact 
by Case A mass transfer if the initial orbital separation is less 
than about 20-30 Rq . They first interact by Case B mass trans¬ 
fer if the initial orbital separation is longer than the boundary 
for Case A mass transfer and shorter than 800-1000 Rq . They 
first interact by Case C mass transfer if the initial orbital sep¬ 
aration is longer than 800-1000 Rq. They do not interact by 
RLOF at all if the initial orbital separation is longer than about 
2000 Rq. The boundaries depend on the initial mass ratio: if 
the mass ratio is larger (for a fixed orbital separation), the pri¬ 
mary overfiows its Roche lobe earlier in its evolution because 
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Fig. 1.— Mass transfer efficiency /3 = AMacc/AMtrans in binaries with 
M\ = IOMq primary stars as a function of initial orbital separation a-yoi nnd 
initial mass ratio ^ini The black, dotted lines indicate the bound¬ 

aries between the three different mass transfer Cases A, B and C. The white 
area hatched from bottom-left to top-right shows binaries that result in stellar 
mergers whereas those hatched from bottom-right to top-left show binaries 
that go through at least one common envelope phase during their evolution. 
Binary evolution is followed until the secondary stars leave the MS. 

the Roche lobe is smaller. The boundaries between the differ¬ 
ent mass transfer cases therefore shift to longer initial orbital 
separations for larger mass ratios. 

There is a small zone/gap above the Case C mass transfer 
region for large initial mass ratios > 0.95) where stars in¬ 
teract by Case C mass transfer, but only after the secondary 
has left the MS, hence the gap. 

In Fig. 1, we also mark those binaries as mergers which 
start RLOF on the ZAMS (we do not treat these binaries in 
our simulations). We assume that stars enter a contact phase 
if their mass ratio at the onset of Case A mass transfer is less 
than 0.56 (Sec. 2.1). Only the products of these Case A merg¬ 
ers are MS stars and only these contribute to our analysis of 
the PDMF. 

There are two more critical mass ratios visible in Fig. 1: we 
assume that stars enter common envelope evolution at the on¬ 
set of Case B RLOF if the donor star is a Hertzsprung gap star 
and the mass ratio is below 0.25. For giant-like donor stars, 
i.e. donor stars with a deep convective envelope (e.g. red su¬ 
pergiants), we use a formula (Eq. 57 in Hurley et al. 2002) to 
calculate the critical mass ratio. Mass loss from a star with 
a deep convective envelope leads to an increase of its radius, 
hence to even more mass transfer and is thus dynamically un¬ 
stable. The giant-like donor star engulfs its companion and 
the binary enters common envelope evolution. If all these cri¬ 
teria do not apply for the binary star at the onset of RLOF we 
use the critical mass ratio = 0.33. 

The transferred mass AMtrans froBi a IOMq primary star 
during RLOF is shown in Fig. 2. The transferred mass is very 
similar in all interacting binaries that do not enter a contact 
phase Case A mass transfer is caused by increasing stellar 
radii due to nuclear evolution during core hydrogen burning. 
MS stars with masses larger than 1.25Mq have radiative en¬ 
velopes in our models. Mass transfer from stars with radiative 
envelopes is stable, i.e. such stars shrink back into their Roche 
lobe as a reaction to mass loss. Case A RLOF starts earlier 
and lasts longer for smaller initial orbital separations, hence 
the smaller the initial orbital separation the larger the mass 
lost by the donor, i.e. the transferred mass AMtrans- During 
Case B and C mass transfer, RLOF is driven by the expansion 
of the stellar envelope. Stars stop overfilling their Roche lobes 
only after losing (nearly) their whole envelope. The envelope 
mass of the IOMq donor stars of the Case B binaries which 
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Fig. 2.— The mass AMtrans transferred from IOMq donor stars by RLOF 
during the MS evolution of the secondary star as a function of the initial mass 
ratio and the initial orbital separation. The hatched regions have the same 
meaning as in Fig. 1. 
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Fig. 3.— As Fig. 2 but for the total mass accreted by the secondary star, 
AMacc, during RLOF. 


do not enter a CE evolution is always about the same, hence 
is the total transferred mass in Fig. 2. 

Figure 3 shows how much of the mass transferred from the 
IOMq donor stars is accreted by the secondary stars during 
their MS evolution. The only limit on mass accretion in our 
models is the thermal timescale of the accretor which ensures 
that the accretors remain in thermal equilibrium, an implicit 
assumption when using fitting functions for single star evolu¬ 
tion. The thermal (Kelvin-Helmholtz) timescale is given by 


TkH 


. ^7 M/MpMenv/MQ 

T^/RqL/Lq ^ ’ 


(9) 


where M is the total mass, Menv the mass of the envelope 
(Menv = M for MS stars), R the radius and L the luminosity of 
the star. Nearly all transferred mass is accreted during Case A 
mass transfer. Mass transfer during the Hertzsprung gap, i.e. 
Case B mass transfer, proceeds on the thermal timescale of 
the primary. At longer initial orbital separations, the thermal 
timescale of the primary is shorter because the primaries have 
larger radii when they overfill their Roche lobes (mass and 
luminosity are nearly constant in the Hertzsprung gap). The 
thermal timescale of the MS secondary at the same time is in¬ 
versely proportional to a power of its mass (tkh M~^ , x > 0). 
We therefore expect less accretion at longer initial orbital sep¬ 
arations (because of shorter thermal timescales of the donor 
stars and hence faster mass transfer) and for smaller initial 
secondary masses (because of longer thermal timescales of 
less massive accretors). This causes the gradient visible in the 
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Fig. 4.— Fraction of simulated binaries that interact by RLOF during the 
MS of the secondary stars and fraction of interacting binaries that merge on 
the MS, that transfer mass stably by RLOF and that go through a common en¬ 
velope (CE) phase as a function of the mass of the primary stars. The number 
of simulated binaries is given by iVsimuiated? the number of interacting bina¬ 
ries by interact 5 the number of binaries that merge on the MS by iVMS merger, 
the number of binaries that transfer mass stably by RLOF by Actable rlof and 
the number of binaries that go through a CE phase by Nce- By interacting 
binaries, we mean that the primary overfills its Roche lobe resulting in mass 
transfer or a merger. 

Case B region of Figs. 1 and 3. 

Combining the results of the transferred and accreted mass 
during RLOF reveals the mass transfer efficiency as shown in 
Fig. 1. Mass transfer is nearly conservative during Case A, i.e. 
almost all transferred mass is accreted (f3 ^ 1), and becomes 
non-conservative during Case B and C mass transfer. 

Next, we analyse the binary parameter space for different 
primary masses. Figures 17-22 in Appendix A contain the 
mass transfer efficiency, transferred and accreted mass as in 
Figs. 1, 2 and 3 but for primary masses of 2, 5, 20, 50, 70 
and IOOMq. Binaries in which the primary star once filled its 
Roche lobe are called interacting binaries. From the data of 
Fig. 1, we compute the ratios of the number of binaries that 
go through a common envelope phase, merge on the MS and 
transfer mass stably by RLOF to the number of interacting 
binaries and the ratio of the number of interacting to simulated 
binaries. We plot these ratios as a function of primary mass 
in Fig. 4 taking into account the initial distribution functions 
of binaries discussed in Sec. 2.2. The overall trend is that the 
more massive the primary star the larger the fraction of stars 
that transfer mass stably by RLOF, the larger the fraction of 
stars that merge on the MS and hence the smaller the fraction 
of binaries that go through a common envelope phase. 

The number of interacting binaries decreases rapidly 
around 50 Mq because more massive stars cross the 
Humphreys-Davidson limit (Humphreys & Davidson 1979) 
after their MS evolution and are subject to strong wind mass 
loss. The mass loss widens the orbits such that binaries can¬ 
not interact by RLOF (cf. Vanbeveren 1991). Consequently, 
the fraction of MS mergers increases because it is normal¬ 
ized by the number of interacting binaries while the number 
of Case B binaries that go through a common envelope phase 
(^ini <0.25) decreases. 

The fraction of interacting binaries shows a kink around 
8Mq because this is the mass above which stars explode as 
supernovae. Our 7Mq model reaches a maximum radius that 
is larger than that of e.g. our SMq model because the latter 
star explodes before reaching a similarly large radius. The 
parameter space for interaction is therefore smaller in bina¬ 
ries with 8Mq primary stars. From thereon, the number of 


interacting binaries gradually increases with primary mass 
because more massive primary stars reach larger maximum 
radii. Above ^ 13 Mq, the maximum radius is larger than that 
of a 7Mq star and so is the fraction of interacting binaries. 

Between 2 and IOMq the fraction of binaries that transfer 
mass stably by RLOF increases while the fraction of systems 
that go through a common envelope phase decreases. Stars 
more massive than about 2Mq expand significantly while 
crossing the Hertzsprung-Russell diagram and the chance of 
interaction by Case B mass transfer is therefore larger. Less 
massive stars are cooler and develop convective envelopes, 
ascending the giant branch (GB), after little expansion. Mass 
transfer from such stars is dynamically unstable leading to 
a common envelope phase. The more massive stars are, the 
larger the Hertzsprung gap, the larger the number of stars that 
interact by Case B mass transfer while having a radiative en¬ 
velope and the fewer the number of binaries that enter a com¬ 
mon envelope phase. 

Stars with M > 13Mq ignite helium in the core during the 
Hertzsprung gap in our models and only slightly climb the 
GB, meaning there is only a limited range of initial orbital 
separations that leads to stars that interact as giants and thus 
enter a common envelope phase. Contrarily, the fraction of 
binaries that transfer mass stably by RLOF plateaus for pri¬ 
mary masses greater than IOMq because the available range 
of separations to transfer mass from a star without a fully con¬ 
vective envelope does not change significantly. 

Binaries in which the primary star is more massive than 
about 22 Mq do not experience Case C mass transfer because 
the primary stars lose their envelopes through strong stellar 
winds during core helium burning which widens the orbit and 
prevents the star from further expansion. As a consequence, 
the star cannot overfill its Roche lobe after core helium burn¬ 
ing and there is no Case C mass transfer. 

Binaries with 0-type primaries, i.e. with Mi > 15Mq are 
distributed according to the initial orbital period distributions 
of Sana et al. (2012). This enhances the number of stellar 
mergers on the MS and the occurrence of Case A mass trans¬ 
fer because there are more binaries with short orbital peri¬ 
ods according to the orbital period distribution of Sana et al. 
(2012) than for a distribution that is fiat in logP as used for 
less massive stars. The change of the period distribution leads 
to the abrupt increase in the fraction of MS mergers around 
15Mq and the decrease of the fraction of binaries transfer¬ 
ring mass stably by RLOF because the emphasis is now on 
binaries in the closest orbits which preferentially merge (see 
Figs. 1 and 17-22). 

We discuss uncertainties in our models regarding single and 
binary star evolution and the initial distribution functions in 
Appendix B. 

2.4. Construction of mass functions 

Constructing the mass function of a population of single 
stars is straightforward: we bin the masses and count the num¬ 
ber of stars per mass bin. In the case of binary stars we con¬ 
struct three different mass functions according to: 

1. The primary masses only. 

2. The secondary masses only. 

3. An observed mass reconstructed from the total stellar 
luminosity. 

In the latter case we derive masses from the total luminosity 
of binaries: we add the luminosities of both components and 
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Fig. 5.— Present-day mass functions of coeval single star populations of 
different ages. Dashed vertical lines correspond to the mass of stars with a 
MS lifetime equal to the age of the population. The most massive stars have 
the shortest MS lifetime and hence disappear first. Stellar wind mass loss 
creates the peaks at the high mass end most prominently seen in the S.OMyr 
old population. 

numerically invert the mass-luminosity (ML) relation given 
by our fitting formula (Eq. 12 of Hurley et al. 2000). We call 
mass functions constructed in this way observed mass func¬ 
tions, because it comes closest to mass functions from pho¬ 
tometry where luminosities are directly translated into masses 
and binaries are unresolved. We emphasize again that we only 
consider MS stars in our analysis. 

3. MODULATION OF MASS FUNCTIONS BY STELLAR EVOLUTION 
3.1. Single star populations 

Two phenomena infiuence the mass functions of a popula¬ 
tion of single stars. First, stars have shorter lives the more 
massive they are. When stars leave the MS, they disappear 
from our mass functions. Second, massive stars lose mass 
through stellar winds. This reduces stellar masses and hence 
changes the mass function. In Fig. 5 we show PDMFs of 
coeval single star populations of different ages. The mass 
functions are constructed from the theoretically known stel¬ 
lar masses (approach 1 in Sec. 2.4). The dashed vertical lines 
mark the initial mass of stars that have a MS lifetime equal to 
the age of the population. Over-plotted in black is the initial 
distribution of stellar masses, i.e. the Salpeter (1955) IMF. 

The most massive stars leave the MS first which causes a 
truncation close to the dashed vertical lines. From here on 
we call this truncation the turn-off because stars at the trun¬ 
cation of the mass function are located at the turn-off in the 
Hertzsprung-Russell (HR) diagram of a star cluster. Stars 
at the turn-off are called “turn-off stars” and their mass the 
“tum-off mass”. Younger than about lOMyr the tum-off and 
the dashed lines do not coincide, they are displaced from 
each other. The displacement is caused by stellar wind mass 
loss and amounts to ~ 40 Mq in the 3.0Myr PDMF for our 
adopted mass loss recipes. 

We use a maximum initial stellar mass of IOOMq. Such 
stars have MS lifetimes of about 3.0Myr, so the upper mass 
boundary plays no role in populations older than 3.0Myr be¬ 
cause all stars with initially masses larger than IOOMq have 
left the MS. Our results are infiuenced by boundary effects at 
ages younger than 3.0Myr. 

A prominent feature in the PDMF is the accumulation of 
stars close to the turn-off because of mass loss through stellar 
winds. The more massive a star the more mass is lost during 


its MS evolution. Stars with masses close to the turn-off mass 
are the currently most massive stars in the cluster and will 
soon leave the MS, so the greatest accumulation of stars in the 
PDMF is found close to the tum-off. In our calculations, stars 
initially less massive than 9 Mq at Z = 0.02 (tms ^ 30Myr) 
lose less than 1% of their initial mass via winds during the 
MS — the accumulation of stars in the PDMF disappears to¬ 
tally and the initial and the present-day tum-off masses are 
the same (the dashed vertical lines coincide with the turn-off). 
Note that we do not see a sharp truncation as indicated by 
the dashed vertical lines because of the finite size of our mass 
bins. 

The magnitude of the wind mass loss peak in the PDMF 
depends on how many stars are shifted to lower masses and 
by how much. It therefore depends on the strength of stellar 
wind mass loss, i.e. also on the metallicity, and the slope of 
the IMF. In Fig. 6 we show this dependence at the high mass 
end of PDMFs (> IOMq) on a linear scale where the size of 
the bump is more apparent: the higher the metallicity, i.e. the 
stronger the wind mass loss, the bigger the peak. Similarly, 
the fiatter the mass function, i.e. the more high mass stars exist 
compared to lower mass stars, the bigger the peak. The num¬ 
ber of stars in the highest-mass bin of the 3 Myr mass function 
is increased compared to the IMF by 127% for F = -1.35 and 
by 162% for F = -0.9 at a metallicity Z = 0.02 and by about 
55% for F = -1.35 and 70% for F = -0.9 at Z = 0.01. 

The wind mass loss peak allows us to determine the overall 
mass lost by stars during their MS evolution. In practice, it 
may be necessary to model luminosity functions rather than 
mass functions with different mass loss recipes to match ob¬ 
served luminosity functions. Only by doing so can one over¬ 
come the inherent problem that masses derived from observed 
luminosities rely implicitly on the mass loss prescription used 
to derive the applied ML relations. 

Figures 5 and 6 show that the wind mass loss peak is only 
visible in sufficiently young clusters in which stellar winds are 
strong. We now investigate the age and metallicity range of 
stellar populations whose mass functions are likely to show 
the wind mass loss peak. Let Mto,p be the present-day tum- 
off mass in the star cluster and N\ the number of excess stars 
in the peak compared to the IMF which is a power-law with 
slope r (Eq. 4). We can then redistribute these stars such 
that the peak is removed by filling up the IMF from the top 
end, thereby obtaining the initial mass of the turn-off stars, 
Mto,i. Let AM be the mass lost by turn-off stars over their 
MS evolution, i.e. AM = Mto,i -Mto,p. Because the number of 
excess stars Ni equals the number of stars N 2 between Mto,p 
(Fig. 6), we have. 


Ni=N2 = N 2 


= -M^ 

r^^^to,p 


( 


Ml 


:) 


r 



( 10 ) 


To judge whether the accumulation of Ni stars is enough to 
see a wind mass loss peak, we have to compare Ni to the 
number of stars expected from the IMF in a mass range that 
is slightly less massive than the turn-off, Aimf. We define Aimf 
to be the number of stars in the mass range from /Mto,p to 
Mto,p where / < 1. For / = 0.75, we require that Ni has to 
exceed 10% of Aimf to be able to see the wind mass loss peak. 
The 10% requirement and the value for / are chosen such that 
a visible wind mass loss peak is predicted in the 7 Myr but 
not 11 Myr stellar populations at Z = 0.01 and Z = 0.02 (cf. 
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Fig. 6.— High mass end (> IOMq) of mass functions of coeval stellar populations with IMF slopes F and metallicities Z. The panels show (a) a population 
drawn from a Salpeter IMF F = -1.35 at solar metallicity Z = 0.02, (b) as (a) with Z = 0.01, (c) as (a) with F = -0.9 and (d) as (c) with Z = 0.01. The black 
dotted line shows the IMF and the filled regions labelled A^i indicate the accumulation of stars because of stellar wind mass loss at an age of 3Myr. Because the 
area labelled Ni contains those stars that depopulated the mass function for masses beyond the mass function peak, it contains the same number of stars as the 
area labelled N 2 . The wind mass loss peaks are stronger the flatter, i.e. the more positive, the IMF slope and the stronger the wind mass loss (i.e. the higher the 
metallicity). 


Fig. 6). Analogously to Eq. (10), we compute Aimf and find 
that, 


Ni 

^imf 


(l+AM/Mto,p)’^-l 

l_/r 


> 10% for/= 0.75 


( 11 ) 


in order for the wind mass loss peak to be visible. This crite¬ 
rion translates into a relative wind mass loss of. 


AM 

^to,p 


A^i 

.Mmf 






- 1 . 


( 12 ) 


According to our criterion, the minimum relative mass loss is 
3.7% for r = -1.35 and 3.3% for E = -0.7. 

In Fig. 7, we show how much mass is lost by stars through 
stellar winds during their MS evolution as a function of initial 
mass and metallicity. We further indicate the corresponding 
MS lifetimes for Z = 0.02 on the top and plot the criterion 
from Eq. 11 for Ni /Mmf = 5,10 and 15%. The mass functions 
of stellar populations of ages and metallicities to the right of 
the A^i/Mmf = 10% line are likely to show the wind mass loss 
peak. As a rule of thumb, the mass functions of stellar pop¬ 
ulations younger than about lOMyr at Z = 0.02 are likely to 
show a visible wind mass loss peak. Note that our criterion 
is a rather rough estimate and depends on the wind mass loss 
prescription and its metallicity scaling. 


3.2. Binary star populations 

Two processes additionally infiuence the mass functions 
of binary star populations compared to those of single stars. 
First, stars can merge during their evolution and, second, 
stars can accrete mass by RLOF. In what follows we discuss 
three different mass functions to disentangle these effects: we 
present mass functions of primary stars, of secondary stars 


Main-sequence lifetime (Z=0.02) / Myr 
24.5 16.6 11.8 8.9 7.0 5.7 4.9 4.3 3.8 3.3 3.0 


0.1 1 


0.01 

1.0 1.2 1.4 1.6 1.8 2.0 

log( Mi.i/Mq) 

Fig. 7.— Fractional initial mass loss, AM/Mini, of stars during the main 
sequence as a function of initial mass and metallicity Z. The solid and dotted 
lines show the criterion of the visibility of the wind mass loss peak in mass 
functions from Eq. (11) for M/iVimf = 5, 10 and 15%, respectively, and an 
IMF slope of F = -1.35. According to our criterion, a wind mass loss peak is 
likely visible in stellar populations to the right of the solid, 10% line (see text 
for more details). 



and those constructed from the ML inversion method (ob¬ 
served mass functions) as discussed in Sec. 2.4. 

The primary star is initially more massive than the sec¬ 
ondary star, so it evolves faster and is the donor star during 
RLOF. Analogously, the secondary star is the mass gainer. In 
the top panel of Fig. 8 we show the mass functions of pri¬ 
mary stars, i.e. stars that have not yet interacted, that have lost 
mass by mass transfer or that have merged. The mass func¬ 
tions in the middle panel of Fig. 8 show secondary stars that 
have not interacted yet and that have gained mass by RLOF. 
The observed mass function in the bottom panel of Fig. 8 is a 
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log M/Mq 

Fig. 8.— Temporal evolution of PDMFs of all primary stars and mergers 
(top panel), all secondary stars (middle panel) and for observed masses for 
all stars (bottom panel; see Sec. 2.4). By construction, merged stars appear 
in the PDMFs in the top panel and mass gainers of RLOF in the PDMFs in 
the middle panel. The vertical dashed lines indicate the mass of stars having 
a MS lifetime equal to the age of the population, i.e. they indicate the initial 
mass of the turn-off stars. Contrary to a population made only of single stars, 
stars more massive than the initial mass of the turn-off stars exist. The black 
solid lines are the Salpeter IMF and the initial distribution function of the 
secondary stars (Eq. 13) respectively. 


combination of all stellar evolutionary effects. 

The primary stars are initially distributed according to the 
Salpeter IMF (Eq. 4) and the secondary stars according to 


din m 2 


iuu p 

/ '0(lnmi)0(lnm2)dlnmi • / x(ln(2)dlna 
Jlnmi J\na 

=1 

(100’"“* -OT2“') , (13) 


where V^(lnmi), 0 (lnm 2 ) and xClna) are the initial distribution 
functions of the primary and secondary masses and the orbital 
separations as defined in Sec. 2.2 (we assume ^min ~ 0.0 in the 
last step). Compared to the IMF of primary stars, the initial 
mass distribution of secondary stars is lowered by a factor 
(1 -r)“^ (?^ 0.43). The decline of the secondary IMF at the 
high mass end is caused by the maximum initial mass (cf. 
the first term on the right-hand side of Eq. 13: the maximum 
initial mass in our calculations is IOOMq). The overall slope 
is F = -1.35 again because of the fiat mass ratio distribution. 

In single stars, no star is more massive than the turn-off. 
Including binary stars it is possible to populate a tail of stars 
that extends the high mass end of single star PDMFs by about 
a factor of 2 in mass. 

In the top panel of Fig. 8 only stellar mergers populate the 
PDMF tail (i.e. the PDMF on the right-hand side of the verti¬ 
cal dashed lines). The relative number of mergers in the tail is 
lower the older the population as expected from our analysis 


of the binary parameter space in which the fraction of binary 
stars that merge on the MS increases with the initial mass of 
the primary star (Fig. 4). 

The maximum mass of stellar mergers is reached by ini¬ 
tially equal-mass binaries in which the final mass is the total 
mass of the binary minus a fraction of 10% that we assume is 
lost during a merger (Sec. 2.1). This corresponds to a mass 
increase by a factor of 1.8 or a shift of about 0.26 dex on the 
logarithmic mass scale. The PDMF at 10Myr extends slightly 
further because of the rejuvenation of the merged stars. The 
mass of the merged star is not increased by more than a fac¬ 
tor of 1.8 relative to the primary mass, but fresh hydrogen is 
mixed into its core. This decreases the fraction of burnt fuel 
and hence the apparent age of the star. Compared to genuine 
single stars of comparable mass, i.e. stars which have not in¬ 
teracted, the stellar merger has more available fuel and thus 
stays longer on the MS. Rejuvenated stars can thus appear to 
be shifted by more than 0.26 dex because the turn-off mass 
decreases simultaneously. 

The PDMF tails in the middle panel of Fig. 8 contain sec¬ 
ondary stars that have accreted mass. The PDMFs again ex¬ 
tend to slightly larger masses than expected from mass accre¬ 
tion alone because of rejuvenation. Mass accretion in young 
stellar populations (ages < 10 Myr) forms PDMF tails that 
even exceed the initial distribution of the secondary stars. At 
later times, the relative number of stars in the tail is less be¬ 
cause the number of interacting binaries that transfer mass 
stably by RLOF decreases for initially less massive primary 
stars (Fig. 4) as does the overall mass transfer efficiency, /3, 
which is coupled to the thermal timescale of the mass gainers 
during Case B mass transfer. The thermal timescales of stars 
become more comparable for larger masses because the ML 
and mass-radius relations are less steep the larger the mass. 

In populations younger than 3.0Myr, the effect of mass ac¬ 
cretion on the PDMF is modest because the most massive 
stars in our models have just left the MS and hence there 
is no contribution from Case B mass transfer to the PDMF 
(Case C does not occur for primary masses larger than about 
22 Mq, see Sec. 2.3). The number of interacting binaries 
drops strongly for primary masses larger than about 50 Mq 
because of the Humphreys-Davidson limit (Sec. 2.3). This 
is important in stellar populations younger than ^ 4.3 Myr 
which is the MS lifetime of 50 Mq stars in our models. 

In the bottom panel of Fig. 8 we use the ML inversion 
method (approach 3 in Sec. 2.4) to construct the PDMFs. 
Both, mergers and accretors show up in the PDMF tails. The 
star with the highest mass in each population is formed by 
RLOF. However, the difference between the highest mass 
achieved by stellar mergers and by RLOF is small — at most 
one bin-width, i.e. Alogm = 0.04 dex — and depends on the 
assumed mass loss in stellar mergers and the mass transfer ef¬ 
ficiency. Only at early times (^3.5 Myr) does the star with the 
highest mass originate in a stellar merger. Unresolved bina¬ 
ries also contribute to the tail of the observed PDMFs but their 
contribution is small (at most 30% of the tail stars are unre¬ 
solved binaries in young, < 3-4 Myr populations) and great¬ 
est for the largest masses because of a flatter ML relation of 
high mass stars compared to lower masses. Unresolved bina¬ 
ries extend the single star mass function by at most ^ 20% in 
mass for 1.. .2M0 stars and by at most ~ 50% in mass for 
60... IOOMq stars (cf. Eq. 15 in Sec. 3.4). 

In Fig. 9 we show the PDMF of secondary stars at three dif¬ 
ferent times — 5, 20 and 100Myr. Case B mass transfer is the 
dominant contribution after 3 Myr, therefore the mass losing 
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Fig. 9.— Mass functions at 5, 20 and 100Myr for secondary stars, i.e. 
mass accretors. The curves show the contributions of the mass transfer 
Cases A, B and C (as defined in Sec. 2.1) to the total PDMF of secondaries 
(black squares). Case B mass transfer is the dominant contribution to the tail 
stars for > 3 Myr, hence the mass loser is not a MS star and therefore not 
present in the PDMFs of primaries in the top panel of Fig. 8. This is why the 
mass gainers are seen clearly in the PDMFs whereas the mass losers are not. 



Fig. 10.— As the bottom panel of Fig. 8, but for a population of stars with a 
primordial binary fraction of /b = 0.5. As discussed in Sec. 3.3 the number of 
stars in the binary tails are halved compared to a pure binary star population. 
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Star has left the MS and is no longer included in the PDMFs 
of the primary stars in the top panel of Fig. 8. The Case A 
mass losers are difficult to find. Nevertheless some of them 
can be found by comparing the PDMFs of populations of sin¬ 
gle (Fig. 5) and binary stars only (top panel of Fig. 8). The 
mass losers are easiest to spot by the difference in the magni¬ 
tude of the accumulation of stars because of wind mass loss 
e.g. in the PDMFs of the 5 Myr old populations. Populations 
made only of single stars show a larger accumulation of stars 
than those made of binary stars, because some donor stars that 
would appear in the wind mass loss peak lose additional mass 
during Case A mass transfer. In practice it is very hard to 
detect the primaries after they have lost mass (de Mink et al. 
2014). 

Case C mass transfer is typically highly non-conservative in 
our models (Sec. 2.3 and Appendix A). The secondary stars 
of such systems therefore do not gain enough mass to make 
any significant contribution to the tail. 

The stars in the tail of the mass functions are rejuvenated 
binary products and hence appear younger than the real age 
of the population — they are blue stragglers. We further char¬ 
acterise these blue straggler stars in Sec. 4 in terms of their 
binary fraction and apparent ages and compare their frequen¬ 
cies to observations. 


log M/Mq 

Fig. 11.— Slopes of the PDMF of single stars (Sec. 3.1) as a function of 
logarithmic mass. Stellar wind mass loss flattens the high mass end of the 
PDMFs. 

(/b = Nb/[Ns +A^b] with Nb the number of binary systems 
and Ns the number of single stars). The considered bina¬ 
ries have initial orbital separations shorter than 10"^ R© — ini¬ 
tially wider binaries are treated as single stars in our models 
(Sec. 2.2). The mass function of a population with this binary 
fraction is then, 

, , 14 ) 

dlogm ydlogm/g ydlogmy^ 

Here (dp/dlogm)^ is the mass function of single stars only 
and correspondingly (dp/d log m)^ is the mass function of bi¬ 
nary stars only. In Fig. 10 we show the observed mass func¬ 
tions of stellar populations with a primordial binary fraction 
of /b = 0.5 (i.e. two out of three stars initially bom in bina¬ 
ries). In comparison to Fig. 8, the number of stars in the bi¬ 
nary tail in the PDMFs (Sec. 3.2) is attenuated. The tail of the 
mass functions of stellar populations with a binary fraction of 
50% is decreased by 0.3 dex, i.e. by a factor of 2, compared to 
those with a binary fraction of 100%. 


3.3. Stellar populations with varying binary fractions 

We use the PDMFs of populations of single and binary 
stars presented in Secs. 3.1 and 3.2 to build stellar popula¬ 
tions composed of a mixture of coeval single and binary stars. 
Fet /b be the binary fraction at birth, i.e. the number of bi¬ 
nary systems divided by the number of total stellar systems 


3.4. Quantification of evolutionary effects on the PDMF 

The stellar wind mass-loss peak at the high mass end flat¬ 
tens the mass function. We quantify the fiattening by com¬ 
puting the PDMF slopes as a function of stellar mass. We 
fit straight lines piecewise to three mass bins at a time by a 
Fevenberg-Marquardt method (Fevenberg 1944; Marquardt 
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Fig. 12.— Ratio of the PDMFs to IMFs for different ages and primordial 
binary fractions. The PDMFs are constructed from the known stellar masses 
(approaches 1 and 2 in Sec. 2.4). 

1963) and show the slopes P of the single star PDMFs in 
Fig. 11. The errors are statistical Icr deviations from the 
best fit. The more massive a star the more mass is lost by 
stellar winds, hence the accumulation and resulting flatten¬ 
ing of the PDMF is strongest in the 3Myr population. The 
slope at the high mass end of the mass function of young 
stellar populations (< lOMyr) is much shallower than the 
Salpeter IMF, F = -1.35. We find extremes of the PDMF 
slope of -0.45 ±0.20 at the highest masses around 5 OM 0 
(logM/M 0 1.7 dex) in the 3Myr population. The flatten¬ 
ing is still significant for lOMyr populations at masses more 
than IOM 0 for which the slope fiattens to values of up to 
P « -0.95 ± 0.09. In other words, not accounting for wind 
mass loss when determining the slope of the IMF for the most 
massive stars may lead to IMF slopes that are biased by up to 
1 dex. 

We further quantify the wind mass loss peak and the bi¬ 
nary tail by dividing the PDMFs of populations consisting 
purely of single stars, purely of binary stars and of a mixture 
of single and binary stars by the initial distribution functions 
(Eqs. 4 and 13) in order to determine the relative importance 
of the evolutionary effects compared to the initial distribution 
of stellar systems. First, we construct PDMFs from the theo¬ 
retically known stellar masses and not using the ML inversion 
method and show the ratio of PDMF to IMF in Fig. 12. This 
enables us to distinguish between the effects of single and bi¬ 
nary star physics and of unresolved binaries on the PDMF. 
Again, the truncations of the PDMFs are due to finite stellar 
lifetimes. 

The PDMF of the 5Myr single star population reaches a 
level of more than 120% of the IMF because of an accumula¬ 
tion of stars caused by wind mass loss (blue solid line close 
to the turn-off mass of logm ^ 1.5 dex). The wind mass loss 
peak becomes weaker with age until it completely disappears 
in stellar populations older than 30Myr (at a metallicity of 
Z = 0.02 with a Salpeter IMF; see also Sec. 3.1). 

Stellar mergers and mass transfer by RLOF shift stars to¬ 
ward higher masses. The number of stars that are slightly 
less massive than the turn-off stars is therefore less than the 
initial number of stars, i.e. PDMF/IMF < 1 (with no wind 
mass loss). The tail of the 5Myr binary population reaches 
an average level of more than 30% of the initial distribu¬ 
tion function, i.e. on average about one third of the IMF in 
a mass range from about 40M 0 (logm ^ 1.6dex) to 80M 0 
(logm 1.9dex) is re-populated by binary evolution. This 


Fig. 13.— As Fig. 12 but constructed using the mass-luminosity inversion 
method (approach 3 in Sec. 2.4). 


level gradually decreases the older the population and the 
smaller the binary fraction. 

Being unable to resolve binaries shifts the observed mass 
function to larger masses (Fig. 13). This is well known (Sagar 
& Richtler 1991; Kroupa et al. 1993; Maiz Apellaniz 2008; 
Weidner et al. 2009) and is evident from the initial mass dis¬ 
tribution of pure binary star populations (OMyr binary popu¬ 
lation in Fig. 13). In a single star population there is no dif¬ 
ference between the observed IMF and the Salpter IMF. With 
unresolved binaries, the PDMF exceeds the IMF and the ra¬ 
tio of PDMF to IMF increases with larger binary fractions. 
In stellar populations with /b =0.5 there are 5 -10% more 
stellar systems than expected from the Salpeter IMF and up 
to ^ 20% more for pure binary populations (/b = 1). The in¬ 
crease is mass dependent: it is greater for larger masses and is 
understood from the ML relation, which can be approximated 
as L = LqM^. In unresolved binaries with equal masses M the 
observed mass Mobs is 

Mobs = ^ = 2*/^, (15) 


i.e. larger than M by a factor of 2^/^. Low mass stars have a 
larger exponent x than high mass stars, hence the larger the 
mass M the larger the factor by which the observed mass is 
increased (e.g. L ex for 1 < M/Mq < 2 and L (x 
for 60 < M/M 0 < 100 ZAMS stars. Tout et al. 1996). This 
translates into an increase of the mass function by a factor of 


V»(lnMobs) 


(16) 


where F is the IMF slope. Unresolved binaries fiatten the 
PDMF most strongly at the largest stellar masses (under the 
assumption that F is constant in the considered mass range). 
We discuss the effect of unresolved binaries on mass functions 
with respect to previous work in more detail in Appendix C. 

Figure 13 shows that unresolved binaries fiatten the PDMF 
at masses less than the turn-off mass instead of steepening the 
PDMF as was the case for resolved binaries. The wind mass- 
loss peaks in the PDMF reach a level of up to ^ 140% of the 
IMF. The tail of the 5 Myr binary population is re-populated 
by on average ± 40% and with maximum levels of about 55% 
of the IMF. The number of stars in the binary tail is halved in 
populations with a binary fraction of /b = 0.5. 

We have shown above that stellar wind mass loss, binary 
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products and unresolved binaries re-shape the high mass end 
of mass functions which may complicate IMF determinations. 
To avoid biased IMF slope determinations in young as well 
as old clusters, we suggest to either exclude those parts of the 
mass function that are expected to be affected by stellar winds, 
binary products and/or unresolved binaries or to correct for 
these effects. As evident from Fig. 13, this concerns roughly 
the mass ranges from 0.5-0.6 to twice the turn-off mass in our 
models (about ±0.3 dex around the turn-off). 

4. BLUE STRAGGLER STARS 

The stars in the tail of mass functions are mainly rejuve¬ 
nated binary products — they are classical blue straggler stars. 
We define a blue straggler as a star whose mass is larger than 
that of the turn-off stars. This definition does not include all 
blue stragglers because there are also rejuvenated binary prod¬ 
ucts that are less massive than the turn-off. In the following 
sections we characterise the blue straggler stars in the binary 
tail and make predictions about their frequencies, binary frac¬ 
tions and ages as a function of cluster age and compare to 
observations to test our predictions and to finally improve our 
understanding of binary evolution. 

Because our stellar evolution code cannot evolve stars more 
massive than IOOMq (Sec. 2.1), boundary effects for ages 
< 3 Myr are expected. This age range is indicated by hatched 
regions in Figs. 14, 15 and 16. 

4.1. Expected and observed blue straggler star frequencies 

Blue straggler stars and their frequencies have been inves¬ 
tigated in the past using population synthesis calculations in¬ 
cluding binary stars (e.g. Collier & Jenkins 1984; Pols & Mar- 
inus 1994; van Bever & Vanbeveren 1998; Hurley et al. 2001, 
2005; Chen & Han 2009; Chatterjee et al. 2013). Qualita¬ 
tively, population synthesis calculations show that binary star 
evolution forms stars that appear as blue stragglers in colour- 
magnitude diagrams through binary mass transfer and merg¬ 
ers. Quantitatively however, the picture is more complicated. 
Some models predict blue straggler frequencies in agreement 
with observations (e.g. Collier & Jenkins 1984; Pols & Mar- 
inus 1994; Hurley et al. 2001, 2005; Chatterjee et al. 2013) 
while others disagree (e.g. Chen & Han 2009; Geller et al. 
2013). This problem is not yet fully resolved. 

Ahumada & Lapasset (2007) present a catalogue of blue 
stragglers in Galactic open star clusters of various ages. They 
count the number of blue straggler stars and the number of 
stars down to two magnitudes below the turn-off based on 
colour-magnitude diagrams. This is a complicated task be¬ 
cause of observational uncertainties such as field contamina¬ 
tion. In addition, the turn-off and age of the cluster are un¬ 
certain, in particular for the younger cluster where only a few 
stars define the turn-off. 

In our models, the number of blue stragglers, Abss, is given 
by the number of stars more massive than the turn-off mass. 
The number of stars two magnitudes below the turn-off, A 2 , 
is computed from the number of stars with luminosities in the 
range to Lto, where Lto is the turn-off luminosity (the 

factor 10“^/"^ corresponds to two magnitudes). 

Sills et al. (2013) show that the number of blue straggler 
stars selected from their position in an observed HR diagram 
following the prescription of Leigh et al. (2011) can be less 
by a factor of about two than the number of blue stragglers 
selected from models based on stellar masses. The reason for 
the difference is that the criteria of Leigh et al. (2011) that 



Fig. 14.— Ratio of the number of blue straggler stars, Abss, to the number 
of stars two magnitudes below the turn-off, A 2 , as a function of cluster age, 
Tci. The raw observational data for each Galactic open cluster of the blue 
straggler catalogue of Ahumada & Lapasset (2007), ALOV, is given by the 
grey triangles, while the same data binned by cluster age is given by blue 
squares. Note that no blue straggler is found in 66% of the clusters younger 
than 500 Myr. These clusters are not on the logarithmic scale but cause the 
decreasing blue straggler frequency of the binned observational data with 
younger cluster ages. Our model predictions are overlayed in black dots. The 
grey regions indicate ages for which boundary effects play a role. 

are used to identify blue stragglers observationally in an HR 
diagram of a modelled star cluster do not cover all blue strag¬ 
glers produced in that model (cf. Fig. 1 of Sills et al. 2013). 
Ahumada & Lapasset (2007) identify blue stragglers differ¬ 
ently in HR diagrams. They count stars as blue stragglers that 
lie between the ZAMS and an isochrone appropriate for the 
cluster. It is therefore likely that the number of observation- 
ally and theoretically identified blue stragglers do not differ 
by a factor of 2 in our case — we expect the difference to be 
smaller but can not rule out systematic differences. Sills et al. 
(2013) further show that there is a strong correlation between 
the number of observationally and theoretically identified blue 
stragglers such that the relative number of blue straggler stars 
as a function of cluster age may be compared directly. 

From the catalogue of Ahumada & Lapasset (2007), we put 
the number of blue straggler stars, Abss, and stars down to 
two magnitudes below the turn-off, A 2 , in age bins of size 
0.25 dex and plot the ratio Abss/Az as a function of cluster age 
in Fig. 14. Additionally, we show the individual data points 
for every cluster and our model predictions. In some clusters, 
no blue straggler star is found. These clusters are not visible 
on the logarithmic scale used in Fig. 14 but do contribute to 
the binned data. 

In our models, the ratio Abss/Az decreases in older popu¬ 
lations. Binary evolution is more efficient in producing blue 
stragglers in high mass, young binaries than in low mass, old 
binaries. This is because the number of binaries that produce 
blue straggler stars by MS mergers and by stable mass trans¬ 
fer, i.e. the number of binaries that do not go through a com¬ 
mon envelope phase, is larger in high mass than in low mass 
binaries (cf. Figs. 4 and 17-22). The peak in the predicted ra¬ 
tio Abss/A 2 around 10 Myr (logTd/yr = 7) is due to a change 
in the distribution functions of the initial orbital periods from 
Opik’s law to the results of Sana et al. (2012) for 0-star bina¬ 
ries (Sec. 2.2). The orbital period distribution of Sana et al. 
(2012) favours tight over wide binaries more than Opik’s law 
does and thus binaries that form massive blue straggler stars. 

A similar decreasing blue straggler ratio Abss/Az in older 
clusters is predicted from binary evolution calculations by 







12 


F.R.N. Schneider et al. 


Chen & Han (2009) but our ratios are larger by factors of 
3-5. Chen & Han (2009) notice a similar difference (of a 
factor of at most two) between their predicted blue straggler 
frequencies and those computed with the population synthe¬ 
sis code of Hurley et al. (2002). Chen & Han (2009) find that 
the difference is related to how they treat binaries in counting 
the number of stars two magnitudes below the turn-off. We 
count each binary as one stellar system because binaries are 
unresolved in the colour-magnitude diagrams from which the 
observed blue straggler frequencies are deduced. We notice 
another difference that might explain the discrepancy: we use 
a maximum initial orbital separation of 10"^ Rq while Chen 
& Han (2009) use 5.75 x 10^R©. Hence, our models con¬ 
tain more interacting binaries (cf. Fig. 4) that can form blue 
straggler stars. In other words, our effective binary fraction is 
larger. 

The observations indicate an opposite trend to our mod¬ 
els: the ratio A^bss /^2 increases with cluster age. There are 
complications that render a direct comparison between the 
observed and our predicted ratio A^bss/A ^2 difficult. There ex¬ 
ist two classes of star clusters in the catalogue of Ahumada 
& Lapasset (2007): clusters with and without blue straggler 
stars. The bi-modality is mainly found in clusters younger 
than about 500 Myr. There is no blue straggler star in 3% of 
their star clusters older than 500Myr, in 38% of star clusters 
with ages between 100 and 500 Myr and in 80% of star clus¬ 
ters younger than 100 Myr. This bi-modality is not understood 
but causes the observed decreasing blue straggler frequency 
with younger cluster ages. 

Stochastic sampling (Schneider et al. 2014b and Ap¬ 
pendix D), i.e. lower number statistics in young clusters com¬ 
pared to older clusters because of the IMF, cannot explain this 
bi-modality. Stochastic sampling is expected to increase the 
scatter between the blue straggler frequencies of clusters of 
similar ages but not to create a bi-modal distribution. How¬ 
ever, stochastic sampling may possibly complicate the accu¬ 
rate determination of the turn-off. 

Supernova kicks (e.g. Lai 2001) might be part of the so¬ 
lution to this discrepancy. The binary companion of a blue 
straggler star that formed by RLOF, i.e. the former mass 
donor, will explode if it is massive enough to undergo core 
collapse. The supernova can break up the binary such that 
the blue straggler leaves the cluster as a runaway star, thereby 
reducing the number of blue stragglers observed in the clus¬ 
ter. Some clusters might lose all their blue stragglers in this 
way, giving rise to the observed bi-modality. The observed bi¬ 
modality is found in star clusters younger than 500 Myr. How¬ 
ever, supernova explosions occur only in young star clusters 
(< 40-50Myr) with massive stars (> 7-8 Mq). Supernova 
kicks can therefore explain only part of the discrepancy. 

The primordial binary fraction in our models is 100%. A 
smaller binary fraction linearly decreases the number of blue 
stragglers and hence the ratio A^bss/A( 2 - Our predictions are 
therefore likely upper limits, implying that our models under¬ 
predict the observed blue straggler frequencies in > 100 Myr 
star clusters and that additional channels for the formation of 
blue stragglers are required. Geller & Mathieu (2011) find 
carbon-oxygen white dwarfs as companions to blue stragglers 
in NGC 188. The white dwarf mass and period distributions 
appear to be consistent with a Case C RLOF formation sce¬ 
nario in which an asymptotic giant branch (AGB) star trans¬ 
fers mass to a MS star. In our models, RLOF from giants typ¬ 
ically leads to common envelope evolution and thus not to the 


formation of blue stragglers. However, Chen & Han (2008) 
show that Case B and C mass transfer from giants to MS stars 
can form blue stragglers. Dynamical cluster evolution also 
produces blue stragglers in stellar collisions (e.g. Geller et al. 
2013) and is expected to be more efficient the higher the den¬ 
sity of the star cluster (e.g. Chatterjee et al. 2013). Another 
form of mass transfer, wind RLOF from AGB stars, might 
be efficient enough to also form blue straggler stars (e.g. Mo- 
hamed & Podsiadlowski 2007; Abate et al. 2013). All three 
contributions, the Case B and C formation scenarios, mergers 
due to collisions and wind RLOF, are missing in our predic¬ 
tions and might potentially help explaining the too low BSS 
predictions in > 100 Myr star clusters (see also Geller et al. 

2013) . 

4.2. Binary fraction of blue straggler stars 

A further testable prediction from our models is the binary 
fraction among blue straggler stars. The blue stragglers that 
form from stable mass transfer by RLOF can potentially be 
observed as binary stars whereas those that form from stellar 
mergers are observed as single stars. It is, however, difficult 
to observationally find the companions of blue stragglers that 
formed by stable mass transfer because the companions can 
be much fainter than the blue straggler (or even be a compact 
object) and the binary orbits are often wide such that most 
searches for radial velocity variations in blue stragglers can¬ 
not detect them (e.g. Collier & Jenkins 1984; Pols & Marinus 
1994). The blue straggler binary might even be disrupted if 
the companion star explodes. In general, post-interaction bi¬ 
naries can often not be identified as such from radial velocity 
variations and mostly appear to be single stars (de Mink et al. 

2014) . 

In our models, the binary fraction of blue straggler stars is 
independent of the primordial binary fraction in star clusters 
because only the sub-sample of primordial binaries can pro¬ 
duce blue stragglers. The binary fraction among blue strag¬ 
gler stars is thus purely determined through binary evolution. 

In Fig. 15 we show the binary fraction /b of blue straggler 
and hence tail stars in our PDMFs as a function of cluster age. 
In our simulations we assume that all binaries with a mass ra¬ 
tio less than 0.56 merge at the onset of RLOF during Case A 
mass transfer (Sec. 2.1), i.e. mass transfer from a MS star. 
Because of the assumed fiat mass ratio distribution, at least 
56% of all Case A binaries merge and are thus observed as 
single stars (/b < 44%). Our analysis of the binary parameter 
space (Figs. 17-22) reveals that more than 56% of all bina¬ 
ries merge during Case A mass transfer because even binaries 
with a mass ratio greater than 0.56 come into contact if they 
are initially in close orbits (e.g. Wellstein et al. 2001). In bina¬ 
ries younger than 2 Myr, only the primary stars of very short 
orbit binaries overfill their Roche lobes. Most of them, re¬ 
gardless of the mass ratio, merge, leading to binary fractions 
< 10%. Later, the binary fraction increases with time because 
also binaries in initially wider orbits interact, producing blue 
straggler stars by stable Case A and B mass transfer without 
merging, i.e. leaving behind a blue straggler in a binary. 

Case A mass transfer predominantly forms blue stragglers 
by merging in more than 56% of all Case A binaries in our 
simulations. Contrarily, Case B mass transfer purely creates 
blue straggler star binaries (neglecting binary disruption by 
supernova explosions). The binary fraction in Fig. 15 there¬ 
fore reaches a maximum of ~ 60-70% and stays at this level 
after all Case B binaries have enough time to interact to con¬ 
tribute blue straggler binaries. 
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Fig. 15.— Binary fraction among blue straggler stars as a function of clus¬ 
ter age, Tci, corresponding to stars in the tail of our PDMFs. The given binary 
fractions are upper limits because we neglect supernova kicks that may dis¬ 
rupt binaries. The grey regions indicate ages for which boundary effects play 
a role. 

In our models, Case B mass transfer in binaries with pri¬ 
mary stars more massive than 5Mq leads to fewer common 
envelope phases and consequently more blue stragglers than 
in binaries with less massive primary stars (cf. Case B regions 
of binaries with 2M0 and 5 Mq primary stars in Figs. 17 
and 18, respectively). The binary fraction of blue straggler 
stars therefore decreases in populations older than the lifetime 
of 5Mq stars, i.e. older than ^ lOOMyr, and is given by the 
fraction of stars that merge during Case A (and early Case B) 
mass transfer. 

Our blue straggler star binary fractions are upper limits be¬ 
cause we neglect supernova kicks that might disrupt binaries 
and because the companions might be hard to detect observa- 
tionally. We exclude supernova kicks because they introduce 
a random process that can only be fully taken into account 
by many repeated calculations which is impractical in our ap¬ 
proach. 

4.3. Apparent ages of blue straggler stars 

Mass gainers and mergers appear to be younger than other 
cluster members because they are rejuvenated by mass accre¬ 
tion (Sec. 2.1). The age of a cluster is given by the MS life¬ 
time of the turn-off stars, rMs(Mo) where x is the 

exponent of the ML relation, L oc M^. The most massive blue 
straggler stars have a mass of about twice the mass of the 
turn-off stars in our models. Their apparent age is thus a fac¬ 
tor rMs(2Mto)/rMs(Mo) = 2^“^ smaller than their true age if 
we neglect mixing of fresh fuel into convective cores which 
makes stars look even younger than the analytic approxima¬ 
tion given here. The ML relation of MS stars is flatter, i.e. 
has a smaller exponent x, at high masses. The exponent ap¬ 
proaches X = 1 in stars close to the Eddington limit and is as 
large as x = 4-5 in low mass stars (^ 1-2 Mq). The most mas¬ 
sive blue straggler with exponent x = 4, i.e. a blue straggler of 
lower mass in older clusters, appears to be rejuvenated by a 
factor 1/8 (the cluster appears eight times as old as the blue 
straggler). Contrarily, the most massive blue straggler with an 
exponent x = 2, i.e. a blue straggler of higher mass in younger 
clusters, appears to be rejuvenated by only a factor of 1/2 (the 
cluster appears twice as old as the blue straggler). 

In Fig. 16 we show the apparent stellar age, of the most 
massive and of the apparently youngest blue straggler star in 
our simulations as a function of cluster age, Td. The above 
mentioned trend for the most massive blue straggler is re¬ 



Fig. 16.— Apparent stellar age, r*, of the most massive and the apparently 
youngest blue straggler star as a function of cluster age, Td. The grey regions 
indicate ages for which boundary effects play a role. 


covered, i.e. the most massive blue straggler in young star 
clusters appears to be less rejuvenated than in old ones: in a 
lOMyr star cluster, the most massive blue straggler appears 
to be younger than the cluster by a factor of 0.48 while this 
factor is smaller, only 0.17, at 1 Gyr. 

The most massive blue straggler is not the one that appears 
youngest (Fig. 16). The apparently youngest blue straggler 
forms from accretion onto a relatively unevolved star which 
then becomes a blue straggler. In our models, the apparently 
youngest blue straggler is about 40% more massive than the 
turn-off and forms from late Case A/early Case B mass trans¬ 
fer in a binary with a mass ratio 0.5-0.6. The mass ratio en¬ 
sures that the secondary star is relatively unevolved. Mass 
accretion then brings the secondary mass above the mass of 
the turn-off, forming a blue straggler close to the ZAMS, i.e. 
apparently very young — the apparent age is of the order of 
7% of the cluster age (cf. Fig. 16). 

The apparent age of the most massive blue straggler de¬ 
pends on how massive blue stragglers can get through merg¬ 
ers and stable RLOF, and on how much fresh fuel is mixed 
into convective cores during these processes. The apparent 
age of the most massive blue straggler is younger the more 
massive the star gets and the more mixing occurs. The appar¬ 
ent age of the seemingly youngest blue straggler depends on 
the choice of the critical mass ratios, ^crit, and the mass trans¬ 
fer efficiency because these parameters determine how small 
the initial mass ratio can be in order to accrete enough mass to 
form a blue straggler star. The smaller the mass ratio, the less 
evolved is the progenitor of the blue straggler and hence the 
younger the blue straggler appears. Note that the mass ratio 
can also not be too small: there has to be enough transferred 
mass such that the secondary mass can exceed the turn-off 
mass. 

5. DETERMINATION OF STAR CLUSTER AGES 

As shown in Sec. 4.3, a stellar population can be 10 times 
older than some of its most massive stars appear to be. Such 
rejuvenated stars can bias the age determination of stellar pop¬ 
ulations because it is not always clear whether a star was influ¬ 
enced by binary mass transfer in the past or not. Especially in 
young clusters with OB stars, the most luminous and hence 
most massive stars are often investigated in detail because 
they are easiest to observe. Their ages are then sometimes 
interpreted as the cluster age although the most massive stars 
are likely rejuvenated binary products. 

To avoid potential biases and confusion of stars with reju- 
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venated binary products such as stellar mergers, stars at least 
half the mass of the most massive cluster members should be 
used to determine ages of stellar populations. This rule of 
thumb is implicitly taken into account when dismissing blue 
stragglers from fitting the turn-off of well-populated star clus¬ 
ters in colour-magnitude or HR diagrams. 

The turn-off may be blurred by binary products and it is 
sometimes, especially in young star clusters, challenging to 
determine its location. In such cases, the mass function may 
provide a promising alternative because possibly rejuvenated 
binary products can be identified by the binary tail and the 
age of the stellar population can be directly determined by 
reading-off the turn-off mass. 

First, we consider single star mass functions to clarify the 
general approach to determine cluster ages from mass func¬ 
tions. Single star mass functions are truncated at the present- 
day turn off mass, Mto,p, because of finite stellar lifetimes. 
The present-day and initial mass of the turn off stars in old 
stellar populations are the same because stars have negligible 
stellar winds. The age of old populations therefore follows 
directly from the MS age of stars with an initial mass of that 
of the turn-off. In young stellar populations, stars lose mass 
by stellar winds and a peak forms in the mass function. The 
present-day mass of the turn off stars that is read-off from the 
truncation of the mass function is no longer equal to their ini¬ 
tial mass because of stellar winds. However, without further 
modeling or the need of stellar models, we can correct for 
wind mass loss by redistributing the number of excess stars 
in the wind mass loss peak, , such that the mass function 
is filled up to the initial mass of the turn off stars, Mto,i (cf. 
Fig. 6 and Sec. 3.1). Rewriting the mass lost by stellar winds 
on the MS, AM = Mto,i -Mto,p, in Eq. (10), we have for the 
number of excess stars in the wind mass loss peak, 

V = ^(Mri-Mrp), (17) 

from which we find the initial mass of the turn-off stars, 

/N p \ 

A^.o,i= , ( 18 ) 

in the cluster. The age of the star cluster follows from the MS 
lifetime of the turn off stars with initial mass Mto,i. 

The mass function of binary stars allows us to determine 
the age of the stellar population in a similar way. The differ¬ 
ence is that the binary star mass functions are not truncated at 
the turn-off mass but rather at about twice this mass. In prac¬ 
tice the mass function may be truncated at less than twice the 
turn-off mass because of stochastic sampling. Consequently, 
another indicator of the turn-off mass than the truncation of 
the mass function is required to determine the age of the stel¬ 
lar population. We use the onset of the binary tail for that 
purpose. In old stellar populations, the onset of the binary 
tail, i.e. the turn-off mass, is indicated by a steep decrease in 
the number of stars. In young stellar populations, the num¬ 
ber of stars in our models do not change that strongly at the 
onset of the binary tail. In such cases, we propose the use of 
the wind mass loss peak to indicate the onset of the binary 
tail and hence to determine the turn-off mass and the age of 
the stellar population. The number of excess stars in the peak 
can be infiuenced by unresolved binaries and binary evolu¬ 
tion (Sec. 3.2). The correction for stellar wind mass loss to 
derive the initial mass of the turn off stars may therefore not 
be possible without the use of stellar models that provide the 


mapping of masses at the end of the MS to initial masses. 

The advantage of this over other methods to determine stel¬ 
lar ages is that it is not biased by apparently younger, rejuve¬ 
nated binary products. Instead, the proposed method identi¬ 
fies and explicitly uses rejuvenated binary products to deter¬ 
mine the turn-off mass and hence to derive the cluster age. 

This new method is used by Schneider et al. (2014b) to de¬ 
termine the ages of the young Arches and Quintuplet star clus¬ 
ters in the Galactic centre. The mass functions of Arches and 
Quintuplet show a wind mass loss peak from which it is possi¬ 
ble to identify the binary tail and derive an unambiguous age. 
The age from the wind mass loss peak results in older cluster 
ages than previously derived for these clusters from the most 
luminous stars. The most luminous stars belong to the binary 
tail of the mass function and are thus likely rejuvenated binary 
products. 

6. CONCLUSIONS 

We use a rapid binary evolution code to investigate how sin¬ 
gle and binary star evolution shape PDMFs with time. To that 
end, we set up coeval populations of single and binary stars, 
follow their evolution in time and construct mass functions at 
ages ranging from Myr up to Gyr. Our code incorporates all 
the relevant single and binary star physics that directly alters 
stellar masses — wind mass loss, mass transfer in binaries by 
RLOF and by winds, stellar mergers and rejuvenation of stars. 

Finite stellar lifetimes truncate the mass functions and wind 
mass loss results in an accumulation of stars in the PDMF 
creating a peak at the high mass end. The magnitude of the 
peak depends on the strength of stellar winds, i.e. on stellar 
mass and metallicity, and the slope of the IMF: the flatter 
the IMF the stronger the peak. The peak can thus be used to 
constrain stellar wind mass loss. We investigate the age and 
mass ranges of stellar populations for which we expect a wind 
mass loss peak in their mass functions. Typically, the peak is 
present in stellar populations younger than about 10 Myr at 
Z = 0.02 corresponding to stars initially more massive than 
about 18Mq (Figs. 6 and 7). Less massive stars have too 
weak stellar winds. The peak fiattens the PDMF slopes at the 
high mass end by up to 60% in 5 Myr old stellar populations 
with a Salpeter IMF (IMF slope F = -1.35; Fig. 11). 

Binary interaction, i.e. mass transfer and stellar mergers, re¬ 
shape the PDMF at the high mass end, forming a tail which 
extends the PDMF of single stars by a factor of 2 in mass. 
The PDMF tail consists of rejuvenated binary products that 
are not expected to exist from single star evolution and are 
better known as blue stragglers. Binary interactions are more 
efficient in producing the binary tail at high masses (Sec. 2.3). 
The number of rejuvenated binary products in young star clus¬ 
ters (^ 5 Myr) reaches more than 30% of the initial number 
of stars in the mass range corresponding to the tail (Figs. 12 
and 13). Binary interactions are therefore efficient in repopu¬ 
lating the high mass end of PDMFs even in star clusters that 
are only a couple of million years old. 

Unresolved binaries flatten the slope of the PDMF. The 
slope is fiattened by about 0.1 in our zero-age populations, 
which is in agreement with previous work on unresolved mul¬ 
tiple systems. Altogether, stellar winds, mass exchange in bi¬ 
nary systems and unresolved binaries reshape the high mass 
end of mass functions within ±0.3 dex of the turn-off mass 
and can therefore bias IMF determinations. 

We compare our predicted blue straggler frequencies to 
those from the blue straggler catalogue of Galactic open star 
clusters (Ahumada & Lapasset 2007). Our models predict 
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a decreasing blue straggler frequencies with increasing clus¬ 
ter age — in contrast to the observations. The observed blue 
straggler frequency drops with younger ages because there are 
no blue stragglers identified in two-thirds of the open clusters 
younger than 500 Myr but only in 3% of clusters older than 
that. Our models predict about the right amount of blue strag¬ 
glers in young clusters (< 100Myr) but too few in older clus¬ 
ters. Additional blue straggler formation channels, such as 
mergers resulting from stellar collisions and wind Roche-lobe 
overflow from AGB stars, are likely required to explain the 
observed frequency of blue stragglers in older clusters. Also, 
the treatment of Case B and C mass transfer may need revi¬ 
sion in our models to allow for the formation of blue stragglers 
from RLOF of giants to MS stars. 

The binary fraction among the blue straggler stars in the 
tail of our PDMFs varies between 40% and 70%, depending 
on cluster age. It is largest in young 10Myr) and smallest 
in old stellar populations 5000Myr). Our binary fractions 
are upper limits because we neglect supernova kicks that may 
disrupt binaries. 

The most massive blue stragglers in the binary tail have ap¬ 
parent ages that are younger by factors of 0.17-0.48 than the 
real cluster age, i.e. cluster ages inferred from the most mas¬ 
sive stars would be too young by factors of 2-6. Some of the 
less massive blue stragglers may show apparent ages that are 
even younger by a factor of 10. Cluster ages derived from the 
most luminous stars should therefore be treated with caution 
because of a likely confusion with rejuvenated binary prod¬ 
ucts. If in doubt, and to avoid potential confusion, cluster 
ages should not be based on the most luminous cluster mem¬ 
bers. Instead, we propose the use of mass functions to identify 
rejuvenated binary products and to derive cluster ages. This is 
possible because the turn-off mass and hence the cluster age 
can be directly determined from mass functions. 

In old stellar populations, the turn-off mass is best indicated 
by a steep decrease of the number of stars at the onset of the 
binary tail. In younger populations, this decrease of the num¬ 
ber of stars is less and the wind mass loss peak may be used 
instead to determine the turn-off mass. Stellar mass functions 
and especially the wind mass loss peak constitute a new and 
unambiguous clock to age-date star clusters. This technique 
is applied by Schneider et al. (2014b) to the Arches and Quin¬ 
tuplet star clusters to determine unambiguous cluster ages by 
identifying likely binary products and thereby resolved the ap¬ 
parent age discrepancies among the most luminous members 
in both clusters. 

The binary products in the tail of young PDMFs have po¬ 
tentially far-reaching consequences. So far, they have mostly 
been neglected when computing the feedback from stellar 
populations. However, the stars in the tail are the most mas¬ 
sive stars in a stellar population and can therefore contribute 
significantly to the ionising radiation, the mechanical feed¬ 
back from stellar winds and supernovae explosions and to the 
chemical enrichment. In the youngest star clusters, binary 
products might even become so massive that they explode as 
pair-instability supernovae, thereby contributing significantly 
to the metal production in the Universe (Heger & Woosley 
2002; Cal-Yam et al. 2009; Danger 2012). 
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APPENDIX 

A. BINARY PARAMETER SPACE CONTINUED 

In Sec. 2.3 we describe how much mass is transferred and 
accreted in our binary models with IOMq primary stars to un¬ 
derstand quantitatively how binary evolution shapes the high 
mass end of PDMFs. Here, we continue this description by 
providing figures equivalent to Figs. 1, 2 and 3 but for pri¬ 
mary masses of 2, 5, 20, 50, 70 and IOOMq (Figs. 18-22, 
respectively). These analyses enable us to fully understand 
the quantitative results presented in this paper. The top panels 
(a) contain the mass transfer efficiency /3 as defined in Eq. (8), 
the middle panels (b) the mass transferred from the primary to 
the secondary stars during stable RLOF and the bottom pan¬ 
els (c) the mass accreted by the secondary stars during stable 
RLOF. 
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Fig. 17. — Mass transfer efficiency (3 (top panel (a)), transferred mass from Fig. 18. — As Fig. 17 but for 5 Mq primary stars, 

primary to secondary during stable RLOF (middle panel (b)) and mass ac¬ 
creted by secondary star during stable RLOF (bottom panel (c)) as functions 
of the initial mass ratio q\ni (i.e. initial secondary mass) and initial orbital 
separation for 2M0 primary stars. The shaded regions have the same 
meaning as in Fig. 1 and indicate binaries which merge and/or go through a 
common envelope phase. 
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Fig. 19.— As Fig. 17 but for 20M© primary stars. 


Fig. 20.— As Fig. 17 but for SOM© primary stars. 
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Fig. 21.— As Fig. 17 but for TOM© primary stars. 


Fig. 22.— As Fig. 17 but for lOOM© primary stars. 
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B. UNCERTAINTIES IN THE MODELS 
B.l. Single star evolution 

Stars more massive than IOMq accumulate in a peak at the 
high mass end of PDMFs because of stellar wind mass loss 
(Sec. 3.1). The magnitude of the accumulation depends on 
the IMF slope F and on the strength of stellar winds, i.e. on 
the wind mass loss prescription and the metallicity. Our wind 
mass loss prescription for MS stars (Nieuwenhuijzen & de 
Jager 1990; Kudritzki et al. 1989) tends to slightly underes¬ 
timate stellar wind mass loss compared to Vink et al. (2000, 
2001). We compare our initial and end-of-MS stellar masses 
to the latest non-rotating stellar models of Brott et al. (2011) 
and Ekstrom et al. (2012) and find that the turn-off masses 
agree within 1-3% for stars less massive than about 50 Mq 
and that we overestimate the turn-off masses by up to 20 % 
for more massive stars. This deviation is primarily because 
of the Wolf-Rayet wind mass loss rates used in Brott et al. 
(2011) and Ekstrom et al. (2012) for massive MS stars which 
we only apply for post-MS stars. Increasing the wind mass 
loss for stars > 50 Mq to match the turn-off masses of the lat¬ 
est detailed stellar models results into stronger peaks at the 
high mass end of PDMF younger than about 4.3 Myr. Older 
populations are not affected. 

The wind mass loss peak fiattens PDMFs (Sec. 3.4). 
Whether this effect needs to be taken into account when de¬ 
riving the IMF from observations depends on how mass func¬ 
tions are constructed from a measured sample of stars. If, on 
the one hand, measured stellar luminosities are converted to 
masses by means of a suitable ME relation, the fiattening of 
the PDMF by wind mass loss needs to be taken into account. 
If, on the other hand, the observed stars are compared indi¬ 
vidually to stellar evolution tracks to find their initial masses, 
this effect does not need to be corrected for because stellar 
tracks usually include wind mass loss. Alternatively, the mass 
function slope can be determined omitting the high mass end 
and therefore the wind mass loss peak. 

B.2. Uncertain binary physics 

In this section we discuss changes in the treatment of bi¬ 
nary interactions and their importance for the PDMFs. Be¬ 
sides mass, angular momentum is transferred during RLOF. 
Packet (1981) found that a uniformly rotating star needs to 
accrete only 5-10% of its initial mass through a disc to reach 
a critical (Keplerian) velocity at the equator such that the out¬ 
ermost layers are no longer bound to the star (see also Petrovic 
et al. 2005). This point of view is debated and there are sev¬ 
eral arguments regarding accretion and decretion disks which 
might be able to dissipate angular momentum such that stars 
stay always below critical rotation and thus can accrete much 
more mass (Lin & Pringle 1976; Popham & Narayan 1991; 
Krticka et al. 2011). If a star is spun up to over-critical ro¬ 
tation it probably sheds as much mass as is needed to rotate 
below critical. Mass is lost from the system and takes away 
angular momentum from the star. In the case of RLOF this 
means that not all transferred mass is accreted, mass transfer 
is non-conservative and typically only a few percent of the 
transferred mass is accreted (Danger 2012). Low mass trans¬ 
fer efficiencies reduce the PDMF binary tail. For stellar merg¬ 
ers we assume that 10 % of the system mass is lost and takes 
away the excess angular momentum. This is likely rather an 
upper limit because detailed collision simulations of massive 
MS stars show that < 10% is lost even in equal-mass merger 
events (Glebbeek et al. 2013). 


Detailed binary models find mass transfer efficiencies of 
about 100% for Case A and of the order of 40-100% for early 
Case B mass transfer (Wellstein et al. 2001). Also, tides in 
short period binaries and magnetic fields in combination with 
mass loss (e.g., by a stellar wind), i.e. magnetic braking, are 
efficient in dissipating angular momentum such that stars do 
not reach critical rotation during RLOF. Magnetic fields seem 
to be even generated by strong shear because of mass accre¬ 
tion (see e.g., Plaskett’s star, Crunhut et al. 2013) and stellar 
mergers (Tout et al. 2008; Ferrario et al. 2009; Danger 2012). 

It is mentioned by Hurley et al. (2002) that the critical 
mass ratio ^crit = m 2 /mi = 0.25 that determines whether stars 
come into contact during thermal timescale mass transfer, 
e.g., when the primary star is a Hertzsprung gap star, is 
rather approximate. Case B mass transfer can lead to a con¬ 
tact phase for extreme mass ratios because the mass trans¬ 
fer timescale becomes faster than the thermal timescale of the 
accretor (e.g. Ulrich & Burger 1976; Kippenhahn & Meyer- 
Hofmeister 1977; Neo et al. 1977; Pols & Marinus 1994; 
Wellstein et al. 2001). We implement a similar criterion with 
the consequence that nearly all binary systems undergoing 
Case B mass transfer go through a contact phase in which both 
stars merge or eject their common envelopes. Merging stars 
results in a post-MS object and hence a reduction of MS sec¬ 
ondary stars. Ejecting envelopes terminates any mass transfer 
such that the initial secondary stars cannot accrete any matter 
— in other words the mass transfer efficiency is 0 % for most 
Case B binaries. The resulting changes to our PDMFs are 
small, because we already limit mass accretion to the thermal 
timescale of the accretors. 

The initial distribution functions of secondary masses (or 
mass ratios) and orbital periods are also important for our re¬ 
sults. We assume an initially fiat mass ratio distribution, i.e. 
all mass ratios are equally probable. Small initial mass ra¬ 
tios often lead to contact phases (see the discussion of the 
binary parameter space in Sec. 2.3). A distribution function 
of initial mass ratios which favours equal mass binaries en¬ 
hances the effect of mass transfer and reduces the effect of 
MS coalescence. Also the Case B contribution, which is the 
dominant mass transfer channel at later times (Fig. 9), is in¬ 
creased because the Case B mass transfer efficiency is higher 
for larger mass ratios. If the initial orbital separation distri¬ 
bution favours small initial orbital separations compared to 
Opik’s law. Case A and B mass transfer (including stellar 
mergers) occur more frequently which increases the magni¬ 
tude of the tail of the PDMFs. 

Different IMF slopes T enhance or reduce the importance 
of binary physics. Binary evolution usually increases the mass 
of stars: let this increase be given by a multiplicative factor a, 
i.e. the increased mass is aM. If we increase all stellar masses 
M by the same factor a, we enhance the mass function at the 
increased mass aM by the factor. 


'^(InM) 
7 /; (In aM) 


(Bl) 


The steeper the slope of the IMF, i.e. the more negative F, the 
bigger the effect of any mass gain. Also the fiattening of the 
PDMF because of unresolved binaries is more important for 
steeper IMF slopes and less important for fiatter IMF slopes 
(Sec. 3.4). This trend is opposite if mass is lost e.g., by stellar 
winds: then, a fiatter IMF enhances and a steeper one reduces 
the accumulation of stars. 
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C. UNRESOLVED BINARY STARS 

Binaries in photometric studies are usually unresolved 
which leads to an overestimation of their masses (see Eq. 15). 
A mass function, even at zero age, thus looks different com¬ 
pared to the IMF. We find changes in the slope of the mass 
function of the order of A/3 « 0.1 toward flattened mass func¬ 
tions for zero age populations. The flattening because of un¬ 
resolved binaries is more important for larger masses because 
of the mass dependence of the power law index of the ML 
relation (Sec. 3.4). The problem of unresolved binaries is 
well known (Sagar & Richtler 1991; Kroupa et al. 1993; Maiz 
Apellaniz 2008; Weidner et al. 2009). To compare with their 
results it is important to know how they distributed stars in 
binaries (see e.g., the review by Bastian et al. 2010). 

Sagar & Richtler (1991), for example, do not use a fiat mass 
ratio distribution, but they draw both binary components ran¬ 
domly from one IMF. For an initial IMF slope of F = -1.5 
(our initial IMF slope is Salpeter, i.e. F = -1.35) and a binary 
fraction of 50% they derive a PDMF slope of ^ = -1.16; for 
a binary fraction of 100% they arrive at /3 = -1.10. These 
changes occur according to their analysis in a mass range of 
2-14 Mq and are larger than we find. The major difference 
between their calculation and ours is the assumed distribution 
function of stars in binaries (mass ratios and IMF slope). 

Maiz Apellaniz (2008) not only investigates unresolved bi¬ 
nary stars but also higher order multiples and chance super¬ 
positions in dense clusters. Taking only the effect of binary 
stars into account they conclude “that for most cases the ex¬ 
istence of unresolved binaries has only a small effect on the 
massive-star IMF slope” of the order of | A/3| 0.2. They use 
the same mass ratio distribution as in our analysis but slightly 
different IMF slopes. 

Weidner et al. (2009) investigate the effect of unresolved 
binaries and higher order multiples for different pairing meth¬ 
ods. None of their pairing methods correspond to our fiat mass 
ratio distribution. A steepening of the observed PDMF by 
A/3 = 0.1 is reported. Their PDMF is steeper because they 
take a certain number of stars given by their IMF and pair 
them randomly into binaries (their RP method). The proba¬ 
bility of finding a massive star grouped with another massive 
star is therefore less than to find it grouped with a lower mass 
companion, hence estimated system masses do not change 
significantly for massive stars. At lower masses it is opposite: 
relatively more binaries with stars of similar masses are found 
and the companions lead to higher system mass estimates. In 
total this leads to a steepening of the PDMF. 

All in all it seems that unresolved binaries have limited in- 
fiuence on the PDMF of zero age populations. Pairing stars 
which are randomly sampled from an IMF steepens the mass 
function, whereas a fiat mass ratio distribution fiattens the 
mass function as is shown in our analysis. Observations of 
massive, i.e. 0-type binaries favour a fiat mass ratio distribu¬ 
tion (Sana et al. 2012), so the high mass end of PDMFs of 
young stellar populations with 0-type stars is expected to be 
fiattened if binaries are unresolved. 

D. STOCHASTIC SAMPLING 

The statement that the most massive stars are likely blue 
stragglers has to be treated with caution, because it is only 
valid in “rich” (i.e. massive) clusters where there are enough 
stars to sample the binary parameter space well — this prob¬ 
lem is known as stochastic sampling. Binary evolution e.g., 
re-populates the 5 Myr PDMF on average by more than 40% 


of the IMF above 40 Mq (logm 1.6). If there is initially 
only one star in this mass range, we expect about 0.4 stars to 
be re-populated by binary evolution — if there are initially 
10 stars (if the cluster is a factor of 10 more massive), we 
expect to find about four binary products. So if the cluster 
is not massive enough, there might be no blue straggler star 
at a certain age. It will also take some time until binary star 
evolution produces the first blue straggler stars. This again 
depends on the cluster richness and on how binaries are dis¬ 
tributed: the shorter the initial period, the earlier the binary 
interaction. When the initially most massive stars in a cluster 
evolve toward the end of core hydrogen burning we know that 
blue stragglers can be produced by MS mergers and Case A 
mass transfer. So the whole contribution of MS mergers and 
Case A mass transfer products of binaries with initial IOOMq 
primary stars to the blue straggler star population is expected 
to be present after about 3 Myr. Blue stragglers can be even 
present at younger ages — this is just a matter of stochas¬ 
tic sampling which is investigated in more detail in Schneider 
et al. (2014b). 
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